Absence of magnetic order in Yb 3 Ga 5 Oi2: relation between phase transition and 
entropy in geometrically frustrated materials 
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From muon spin relaxation spectroscopy experiments, we show that the sharp peak (A type 
anomaly) detected by specific heat measurements at 54 mK for the ytterbium gallium garnet com- 
pound, Yb3GasOi2, does not correspond to the onset of a magnetic phase transition, but to a 
pronounced building up of dynamical magnetic pair-correlations. Beside the A anomaly, a broad 
hump is observed at higher temperature in the specific heat of this garnet and other geometrically 
frustrated compounds. Comparing with other frustrated magnetic systems we infer that a ground 
state with long-range order is reached only when at least 1/4 — 1/3 of the magnetic entropy is 
released at the A transition. 

PACS numbers: 75.50.-y, 75.40.-s, 76.75.+i 



A magnetic phase transition to a long-range order oc- 
curs at low temperature for most crystallographically 
ordered compounds containing a three dimensional lat- 
tice of magnetic ions. Its classical signature is a sharp 
anomaly in the magnetic specific heat (SH) at tempera- 
ture T\ which corresponds to the phase transition tem- 
perature; see e.g. Ref. Q. In this letter, for convenience, 
we label such anomaly a A-anomaly although it may not 
have all the properties attributed to such an anomaly . 

However, for some particular lattice structures, the 
long-range magnetic ordering of the magnetic ions may 
be impeded by their geometric arrangement which gives 
rise to frustration of their magnetic interactions 0- 
Most of the experimental studies focus nowadays on py- 
rochlore and gallium garnet compounds, R2T2O7 and 
R3GasOi2 respectively. R denotes a rare earth atom 
and T a transition element. The R ions are arranged on 
a motif of corner sharing tetrahedra for the pyrochlore 
structure. In the garnet case, the R atoms form two in- 
terpenetrating, non-coplanar, corner sharing triangular 
sublattices. For an experimental review, see e.g. Ref. 

The results from SH measurements are particularly in- 
triguing. Although |#cw|/7a > 1, where 9cw is the 
Curie- Weiss temperature, long-range order may not be 
present below T\. 

For example, let us review data on pyrochlore com- 
pounds. According to SH, Gd2Ti2 07 has two magnetic 
phase transitions at 1 K and ~ 0.75 K and no defined 
structure in the SH data is observed above 1 K |. 
The detection of magnetic Bragg reflections from neu- 
tron scattering at 50 mK shows that a magnetic struc- 
ture is established @. Similarly, Gd2Sn 2 07 has a large 
A-anomaly in SH and Mossbauer spectroscopy @ and 
muon spin relaxation (/^SR) |2( are consistent with the 
presence of long-range ordering. A A-anomaly is detected 



at ~ 1.2 K for Er2Ti207 |lfl an d long-range order is 
observed at low temperature |ll|. For Yb2Ti 2 07 a A- 
anomaly is found at T\ ~ 0.21 K 01 but there are no 
long-range magnetic correlations below T\ but rather dy- 
namical hysteretic short-range correlations 0] . We note 
that beside the A peak, a broad peak centered near 2 K is 
present in SH. In the popular spin-ice systems Ho2Ti207 
and Dy2Ti2 07, only a broad anomaly is present in SH 
|l3| and at low temperature the magnetic moment are 
frozen with no long-range order 

Concerning the garnets, only a broad SH pea k is de- 
tected at low temperature for Gd 3 Ga 5 0i 2 njj. This 
compound does not display any long-range order |l7| . 
Dy 3 Ga50i2 displays an SH anomaly at T\ ~ 0.37 K, 
overlapping with a broad peak 0|. Neutron scattering 
results show the presence of a long-range magnetic or- 
der below T\ [lflj . Finally, while a A-anomaly is found 
for Yb3GasOi2 at T\ = 54 mK ji^J], no information is 
available in the literature about the existence of a short 
or long-range magnetic order. We note that a relatively 
large SH hump is also present, centered at ~ 0.2 K. 

Since no A-anomaly is present in the spin-ice systems 
and in GdaGasO]^ we shall not consider these systems 
any longer. To gain further understanding of the rela- 
tionship between entropy and magnetic correlations in 
Yb3GasOi2, we performed /xSR measurements on this 
system. Here we assume, that the measured SH corre- 
sponds to magnetic dipole moment degrees of freedom 
and not to another exotic order parameter |2lj|. 

In the garnet lattice (space group laid), the crystal 
field acts on the 2 F7/ 2 state of a Yb 3+ ion to leave a 
well isolated ground state Kramers doublet and three 
closely grouped excited doublets having an average en- 
ergy corresponding to A avc /fcB — 850 K .22]. In a good 
approximation, the ground state of Yb 3+ ions diluted in 
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FIG. 1: Zero-field /iSR spectra recorded for Yb3GasOi2, both 
above and below T\. The solid lines in the /^SR spectra are the 
results of fits to exponential relaxation functions as explained 
in the text. In the inset 170 Yb Mossbauer reproduced from 
Ref . [25| . These data show the absence of long or short-range 
correlations below T\. 



Y 3 GasOi2 is described by an effective spin S' = 1/2 with 
an isotropic ^-factor, g = 3.43 jU 0|- The same de- 
scription of Yb 3+ in Yb3GasOi2 is expected to be valid. 
This is supported by the electronic entropy variation be- 
low 1 K which is very close to R In 2 [2(| , as expected for 
a doublet. R is the gas constant. From the measured 
g-factor and bulk magnetization |20l| . the mean value of 
the Yb 3+ magnetic moment at low temperature is found 
to be ~ 1.7 fiB- The interionic interaction is essentially 
antifcrromagnetic as shown from the susceptibility above 
T\ which is lower than that extrapolated from high tem- 
perature for a Curie- Weiss behaviour [ioj . From the po- 
sition of the maximum of the bump of the SH an estimate 
for this interaction strength corresponds to at least 0.2 K. 
The antiferromagnetic nature of the interactions together 
with the crystallographic structure lead to frustration. 

A sample was prepared by heating the constituent ox- 
ides to 1100° C four times with intermediate grindings. 
Its quality was checked by x-ray diffraction and suscepti- 
bility measurements. Zero-field /iSR measurements were 
performed at the ISIS and PSI muon facilities covering 
the temperature range from 21 mK to 290 K. Additional 
spectra were recorded with a longitudinal field. 

170 Yb Mossbauer spectroscopy experiments |25| per- 
formed between 36 mK and 4.2 K show no resolved hy- 
pcrfinc structure at any temperature including below T\ 
as depicted in the inset of Fig. ^ This means that the 
moments fluctuate at a frequency larger than ~ 300 /xs _1 
down to 36 mK. This also indicates the absence of mag- 
netic correlations either short- or long-range. 

We now present zero-field /uSR results. The tech- 
nique consists of implanting polarized (along direction 



Z) muons into a specimen and monitoring P^ xp (£) which 
is the evolution of the muon ensemble polarization pro- 
jected onto direction Z [ifjj . The quantity actually 
measured is the so-called asymmetry corresponding to 
a P| xp (i), where a ^ 0.24 

For a paramagnet, P^ xp (t) tracks the dynamics of the 
magnetic field at the muon site, Bi oc , reflecting the dy- 
namics of the electronic moments. In the fast fluctuation 
or motional narrowing limit, P^ p (t) takes an exponential 
form characterized by a relaxation rate, Xz- A stochastic 
approach to relaxation in zero-field leads to the relation 
Xz = 27^A| f t c , assuming a single correlation time, r c , 
for the magnetic moment dynamics. Azf is the standard 
deviation of the components of Bi oc and 7^ the muon 
gyromagnetic ratio (7^ = 851.615 Mrad s^ 1 T _1 ). The 
motional narrowing limit is valid if j^Azf i~ c <C 1. 

Two examples of /iSR spectra are shown in Fig. ^ All 
the spectra were satisfactorily fitted to an exponential 
function. The temperature dependence of Xz is shown 
in Fig. Since there is no qualitative change in P^ xp (t) 
and since Xz increases monotonously as the temperature 
is lowered, there is neither long-range nor short-range 
magnetic order in YbaGasO]^. This means that the spe- 
cific heat anomaly at T\ does not correspond to the onset 
of a conventional phase transition, in agreement with the 
170 Yb Mossbauer results. This is our first result. 

The fact that we observe an exponential relaxation 
function below T\ is another important point. It im- 
plies that we are in the fast fluctuation limit (we shall 
check this quantitatively below) and therefore the Yb 3+ 
moments continue to fluctuate rapidly down to the lowest 
temperature investigated. This behaviour is very differ- 
ent from that in Yb2Ti207 where the moments abruptly 
slow down below T\ where they are quasi-static [l^ . 

We now comment and interpret the shape of Xz{T). 

As for nuclear magnetic resonance |2jfl , Xz can be ex- 
pressed in terms of the static wavevector dependent sus- 
ceptibility, x(Q)> an( i hnewidth of the quasielastic peak 
of the imaginary component of the dynamical susceptibil- 
ity, r(Q). These two functions are assumed to be scalar, 
consistent with the quasi- isotropic nature of the Kramers 
ground-state doublet. Following Ref. we write 



Xz = 



Mo 7^ 
16tt 2 



fc B T / C(Q) 



X(Q) d 3 Q 
r(Q) (2rr) 3 



(1) 



/io is the permeability of free space and /cb the Boltzmann 
constant. The term C(Q) accounts for the interaction 
of the muon spin with the Yb 3+ effective spins. The 
integration is over the first Brillouin zone. 

For a Heisenberg magnet at high temperature, i.e. if 
the thermal energy is larger than the exchange energy, 
the Q dependence of x(Q) becomes negligible (see e.g. 
Ref. |29j) and it is simply given by the Curie law. Since 
T(Q) is temperature independent in the same limit |29l |. 
we deduce that Xz should be temperature independent. 
This is effectively observed for 0.4 < T < 80 K. 
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FIG. 2: Zero-field muon spin-lattice relaxation rate, Xz, ver- 
sus temperature measured for YbsGasO^. The solid line is 
the result of a fit to a model explained in the main text. 
The two straight dashed lines for T < 0.4 K down to 21 mK 
are guides to the eyes. The specific heat (from Ref. |20j) of 
Yb3Ga5 0i2 is also reproduced. A marked change of slope in 
X Z {T) occurs at T x = 54 mK. 



However, above 80 K, \z(T) decreases steadily as the 
sample is heated. Such a behaviour is sometimes asso- 
ciated with muon diffusion, but it would be quite un- 
usual in an insulating oxide. In fact the decrease of Xz 
is due to the relaxation of the Yb 3+ magnetic moments 
resulting from an Orbach process, i.e. a two-phonon real 
process with an excited crystal-field level as intermedi- 
ate [3(3, HH- A global fit is achieved with the formula 
A^ 1 = A + B me exp[—A e / (fc B T)]. A e is the energy of 
the excited crystal-field level involved. A ~ A^ 1 when 
Xz saturates at low temperature (but still T > 0.4 K). 
The constant B mc models the magneto-elastic coupling 
of the Yb 3+ spin with the phonon bath. The result of 
the fit for 0.4 < T < 290 K is shown in Fig. Taking 
A e /fc B = A avo /fc B = 850 K, we get A = 0.69 (2) fjs and 
B nlc = 250 (70) /is. Combining the value of Xz in the 
range 0.4 < T < 80 K and the value r c ~ 38 ps obtained 
from perturbed angular correlation and Mossbauer data 
in the same temperature interval |25| . we compute A^p — 
0.16 T which is of the expected magnitude. For compar- 
ison A ZF ~ 0.08 T was found for Yb 2 Ti 2 7 above T\ 

M 

Information relevant to the influence of the frustrated 
nature of the magnetic interactions in Yb3GasOi 2 is ob- 
tained from Xz(T) at low temperatures. As shown in 
Fig. |3 A z starts to deviate at ~ 0.4 K from the be- 
haviour expected for the Heisenberg Hamiltonian in the 
high-temperature limit. As we cool down the sample, 
we first observe a mild linear increase of Xz, with slope 
—2.6 (3) fis~ 1 K~ 1 , followed by a sharp increase below T\ 
with slope -21 (7) /zs -1 K _1 . 

In the temperature range T\ < T < 0.4 K, we note that 



the broad SH peak centered at ~ 0.2 K covers the tem- 
perature regime where Xz(T) varies slowly, suggesting 
that both features reflect the same physics. Since there 
are no excited crystal-field energy levels at low energy, 
the hump can not be of a Schottky-type. This extra SH 
reflects short-range correlations among groups of spins, 
see e.g. Ref. Q. This means that the spin dynamics is 
wavevector dependent below ~ 0.4 K. This dependence 
is expected to generate a mild increase of Xz as shown 
by the explicit computation of Paja and coworkers [33j | 
for a simple model. 

Now we turn our attention to the region T <T\. Be- 
ing in an interstitial site, the muon spin can be strongly 
influenced by magnetic pair-correlations. In contrast, 
since a 170 Yb Mossbauer nucleus is embedded in a Yb 
atom which is magnetic, it is mainly sensitive to the 
self correlations. Taking into account that the dynamics 
measured by Mossbauer spectroscopy is not wildly dif- 
ferent above and below T\ j^, we infer that the sharp 
increase of Xz occuring right below T\ is the signature 
of the building up of magnetic pair-correlations. Down 
to the lowest temperature, the /xSR spectra have been 
recorded in the motional narrowing limit, since we com- 
pute 7 m Azf t c ~ 0.02 <C 1 at low temperature, taking 
as an estimate for r c the hyperfine field correlation time 
from Mossbauer (~ 0.3 ns) and computing Azf — 0.08 T 
from the relation Xz = 27^A| f t c . Therefore, the expo- 
nential form of the relaxation is fully justified. 

An experimental estimate for r c can in principle be 
obtained from the analysis of the field dependence of Xz 
assuming that the properties of the system are not mod- 
ified by an external field. In the case of Yb3GasOi 2 we 
have found that a field as low as 0.3 T has a strong in- 
fluence, increasing rather than quenching Xz- We con- 
clude that the field influences the system: this is not 
completely astonishing recalling that a field of 0.7 T in- 
duces a phase transition for the isomorphous compound 
Gd 3 Ga 5 12 0. 

Yb3GasOi2 is therefore a second system with frus- 
trated magnetic interactions, after Yb 2 Ti 2 7 , where no 
long-range order is found below T\. In the pyrochlore a 
sharp first-order transition appears in the fluctuation rate 
of the magnetic moments and below T\ they continue to 
fluctuate slowly (in the megahertz range) at a tempera- 
ture independent frequency. In the garnet we have a new 
scenario: the fluctuations are still rapid below T\ and 
no abrupt change in their frequency is observed. Only 
dynamical short-range correlations build up. The con- 
tinuous rise of Xz is magnetic in origin, and therefore 
excludes a dimer phase to be formed below T\. 

Among the different frustrated pyrochlores and garnets 
with a A peak in SH that we have considered at the begin- 
ning of this letter, i.e. Gd 2 Ti 2 7 , Gd 2 Sn 2 7 , Er 2 Ti 2 7 , 
Yb 2 Ti 2 7 , Dy 3 Ga 5 0i 2 and Yb3Ga50i 2 , some of them 
order and others do not. One can tentatively find a condi- 
tion for a magnetic order to be present. For all these com- 
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pounds the variation of magnetic entropy at low temper- 
ature, say below ~ 10 K, is close to the expected value de- 
duced from the number of electronic degrees of freedom. 
This is R In 8 for the Gd based compounds and R In 2 for 
the others. Inspecting Refs. H H ES E H E S3 , 
the fraction of the total entropy frozen at T\ is roughly 
60 % for Er 2 Ti 2 7 and Dy 3 Ga 5 0i 2 , 40 % for Gd 2 Sn 2 7 
and 35 % for Gd 2 Ti 2 07 whereas it is only ~ 20 % for 
Yb 2 Ti 2 C>7 and Yb3Ga50i 2 . Therefore we estimate that 
a long-range magnetic order is only present when at least 
1/4 — 1/3 of the magnetic entropy is released at T\. 

While the entropy change at the A-anomaly tem- 
perature is about the same for both YbaGa50i 2 and 
Yb 2 Ti 2 07, the Kramers ground-state doublet is approx- 
imately isotropic for the former and strongly anisotropic 
for the latter. This suggests that a high magnetic 
anisotropy impedes the system from exploring different 
configurations, leading to an abrupt slowing down of the 
correlations. Azf is reduced by a factor ~ 14 [T3 . l32| 
when crossing T\ from above for Yb 2 Ti 2 0y whereas this 
reduction is much less (< 2) for YbsGasO^. This is 
again consistent with the anisotropy difference between 
both systems: we expect Azf to be small if the magnetic 
moments are confined to specific orientations. 

In conclusion, according to wisdom, a A-type anomaly 
is indicative of a second order phase transition at T\ for a 
three-dimensional system. However, due to the presence 
of frustrated interactions, only dynamical short-range 
correlations build up below T\ in Yb 3 GasOi 2 . By com- 
parison with results from other three-dimensional frus- 
trated systems, we infer that the long-range order of the 
order parameter does not occur if a too large fraction of 
the entropy is released at a broad bump located above 
T\ . The anisotropy of the system favors a freezing of the 
magnetic correlations. 
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